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Abstract 

The long-lived stau scenario is an interesting possibility at the LHC. We 
study squark mass measurements in this scenario; in particular, we show that 
left- and right-handed squark masses are both measurable. In SUSY events, 
multiple jets are expected, which become a source of combinatorial backgrounds. 
In order to reduce such backgrounds a hemisphere analysis is applied, and we 
discuss mass measurements of squarks in decay modes of qL into qW° or q'W^ 
and also q~R into qB. 



The CERN LHC experiment is now operating and seeking for the physics beyond 
the standard model. The promising candidate of the new physics is the low-energy 
supersymmetry (SUSY), and many studies are devoted for the discovery or the verifi- 
cation of the minimal supersymmetric standard model (MSSM) at the LHC. However, 
such studies should strongly depend on the details of the MSSM. For example, if some 
charged superparticle has a long lifetime enough to escape from the detector, the 
event topology of SUSY signals becomes quite different from the one in the conven- 
tional studies assuming the existence of stable invisible particle. 

We pay particular attention to the case with the long-lived stau. The lightest 
superparticle (LSP) is stable under the .R-parity conservation and thus important for 
phenomenology. For collider study, however, it is more important what the LSP in 
the MSSM sector (MSSM-LSP) is. This is because the MSSM-LSP would not decay 
into the LSP inside the detector if the lifetime of the MSSM-LSP is long enough. 
In a variety of the SUSY breaking models the lightest stau is predicted to be the 
MSSM-LSP. Even though such a charged particle should not be the real LSP, we 
could consider the scenario with other neutral particle as the LSP. One of the viable 
LSP candidates is gravitino, and in such a case the MSSM-LSP, f\, possibly has a long 
lifetimecj Indeed, the typical decay length of f\ becomes longer than the detector 
size, ~ 10 m, if the gravitino is heavier than 0(1) keV. In the following, we consider 
the case that we can treat f\ as a stable particle in collider study. 

We can observe the stau track when its lifetime is enough long. This is the most 
striking signature of the SUSY event in the long-lived stau scenario, and it is useful 
for the discovery of the SUSY at the LHC [2J. The main background of such an 
exotic charged track is muon, however, the distinction of them will be possible by 
using velocity and momentum information on a charged particle [3], 01 [6] . These 
information can be used for measuring the stau mass: In Ref . [5J E] , it is shown that 
the stau mass can be measured with a good accuracy ~ 0(100) MeV. Stau track 
information is also useful for the determination of other superparticle properties and 
the verification of some interesting scenarios [HI El El [9l UHl LTH [121 [13l HH [151 HSl 
[13 HH1 HH1 [201 EH [221 [231 [2H [25] because we can measure the stau's charge, velocity 
and momentum on event-by-event basis. 

An important fact in these studies is that SUSY events in the long-lived stau 
scenario are thought to be background- free if a relevant velocity-cut is applied to f\. 



A scenario with the long-lived charged particle would conflict the cosmology. For example, such 
a particle may spoil the success of the standard big-bang nucleosynthesis. However, the problem can 
be evaded as long as there exists a neutral weekly interacting particle, e.g., gravitino pQ. 
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However, even in such a case, properties of some superparticle are hardly studied. 
One of examples is the left-handed squark mass measurement. A difficulty comes 
from the combinatorial background of the multiple jets. In addition, qi has two 
decay modes into qW° and q'W^ 1 if they are kinematically possible, therefore the 
number of the event of each process is small. These problems make the measurement 
difficult. Although the latter fact is problematic for the mass measurement, once 
we have measured q~i mass in both the decay modes and confirmed two masses are 
identical, this is an important evidence that the squark is really the left-handed one. 
In this sense, we can determine the squark chirality through the mass measurements. 

In this letter, we propose methods for measuring the squark mass. We study 
two sequential decay chains of the left-handed squark, — > qW° — > qrfi and 
qi — > q'W ± — > q'v T T\. In the measurements, it is important to reduce combinatorial 
backgrounds. For this purpose we use a hemisphere analysis, q^ decay into fi is associ- 
ated with a neutrino, z/ r , which is produced from r decay or decay. The momentum 
reconstruction of this u T is another key-point and will be discussed. We also study the 
right-handed squark mass measurement in the decay chain, q R — > qB — > qrf\. 

Now we start our discussion on the squark mass measurement. It is expected 
that a large number of squarks will be produced at the LHC Particularly, the first 
generation squarks, i.e., u and d, would have large production cross sections since 
their superpartners are valence quarks which have large parton luminosities. (We 
employ q as either u or d.) A produced squark will decay to the LSP under the 
i?-parity conservation. The decay chain depends on the MSSM mass spectrum: We 
study the case that the MSSM-LSP is the lightest stau and the mass spectrum in the 
MSSM-sector is 

m q R > m B > m fi, (!) 

mq L > m w± ~ m W o > m fl . (2) 

Such a mass spectrum is indeed realized in the minimal gauge- mediation model [26],|2"T]. 
In the gauge-mediation model the LSP is generally gravitino, and, as mentioned above, 
the stau lifetime is possibly so long. We assume that the decay length of fi is much 
longer than the size of detectors. For the mass spectrum in Eq. ([T]), q R decays into fx 
through the SUSY decay chain 

q R -> qB -> qrf x . (3) 

This decay chain has a dominant branching fraction for ur and d R in general. Although 
another decay mode, e.g., q R — > qW°, may be kinematically possible, such a decay 
only occurs through a neutralino mixing, which is predicted to be small in a variety 
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of SUSY breaking models, and hence would have only a small branching fraction, qi 
decay process is similar to the qn decay, but, for the mass spectrum in Eq. (jSJ), has 
two decay modes because it is charged under the SU{2)l\ 

q L ->■ qW° -> qrr x , (4) 

qi ->■ g'H /± -> gV r fi. (5) 

The branching fraction of q~L — > q'W ± is approximately twice as large as that of 
q L ->■ <jW°. 

Measurements of masses of 5, and W , which are needed in the following 
arguments, have been discussed in earlier works. Several methods are proposed for 
the neutralino mass measurements. For instance, when a SUSY event contains two 
B and both of them decay into the tt\, followed by r decay, v T momenta can be 
reconstructed from the observed missing energy with the approximation that r and 
v T momenta are parallel, and then is determined by peak analysis j6]. Another 
way is to determine as the endpoint of the invariant mass distribution of r-jet (j T ) 
and f\ in the decay process B — > j T v T f\ [20]. In the paper, it is shown that not only 
rug but also m^o can be measured by endpoint analysis. W ± mass measurement has 
been discussed in Ref. [21], and is determined as the endpoint of the transverse 
mass distribution of (u T , fi) system in the decay W ± — » v T f\. Based on these studies, 
we assume in the following that rug, m^o and m^± have been determined. 

Squark masses can be measured by reconstructing their sequential decay processes, 
Eq. (]3]) - Eq. flS}. Let us first review a method for measuring the squark mass [20] . 
To begin with, consider the decay of squark to the lighter stau through the on-shell 
neutralino, x°> followed by r decay, 

Q -> QX° ~> Q Tf i- (6) 

We consider only the hadronic decay mode of the r-lepton, thus, the signal consists 
of one jet (J), one r-jet and the stau. An undetectable v T also accompanies to this 
signal. Nevertheless, we can reconstruct the full kinematics in the process since the 
reconstruction of the 4-momentum of r is possible as follows: When r produced in 
the x° decay is highly boosted, we can approximate that 4-momentum of r and that 
of r-jet are collinear, that is, they can be written in the form 

Pr=rp jT} (7) 

where r is a rescaling factor which should be in the range r > 1. The value of r is 
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determined, requiring that a pair of r and comes from the x° decay, as 

m\ — m 2 - 

r = ^. (8) 

2 Pj T ■ Pfx 

We know now the kinematics in the squark decay of Eq. ([6]), completely. Then, by 
studying the distribution of the invariant mass, 



M q = y/( Pj + rp jT + Pfl f, (9) 

we can determine the squark mass. Information of the mass is imprinted on the peak 
in the distribution. 

The above analysis can be applied to (Jr decay process in Eq. (J3J). By using rug, 
the rescaling factor r is determined on event-by-event basis, assuming that t^ 1 and 
are produced in B decay, as 

2 2 

m% — mf, 

r = r n = ^ 10) 

In early work [20], it has been shown that can been determined by this method. 

One may think that the analysis is applicable to the <jl mass measurement through 
the process in Eq. (TjJ, replacing with r^o, where 



w° 



Zpj T ■ Pn 



(11) 



Though the event topology is quite similar between two processes, Eq. (JBJ) and Eq. (@J, 
there are some difficulties for the study of qt- SUSY event at the LHC would be char- 
acterized by multiple jets which originate from, e.g., squark decay or gluino decay. 
These jets are really burdensome in the squark mass measurement because they be- 
come a source of combinatorial backgrounds. There is no a priori way to know which 
jet is produced from the decay of the squark that we are interested in to measure 
its mass. To make matters worse, the decay — > q\V° may have a small branching 
fraction since there exists another decay mode into q'W ± . Indeed this is true, e.g., 
in the minimal gauge-mediation model, and the branching fraction of q~L qW° is 
roughly 30%. So, comparing to the right-handed squark case, the expected number 
of the signal is small. For finding out such a signal we need to suppress backgrounds 
which mainly come from a combination of jets and also r-jets. 

In order to reduce the combinatorial backgrounds, there is a useful method called 
hemisphere analysis [2Sj - This analysis separates reconstructed objects, such as jets, 
leptons and staus, into two groups. When the analysis successfully works, each group 
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has constituents which arise from the decay of a heavy superparticle produced in hard- 
collision. At first, we need to set a seed for each of two hemispheres. In the long-lived 
stau scenario each SUSY event contains two staus, and we assign one of staus (fi) to 
the seed of the first hemisphere (Hi) and another stau (72) to the seed of the second 
hemisphere (H 2 ). This assignment is quite reasonable because most of stau in SUSY 
event is produced from the decay of a heavy superparticle. Next, the other objects, 
labeled as i (= 1, 2, 3, ...), are assigned to either of new hemispheres, H[ or H 2 , by the 
criterion: 

if d(p Hl ,Pi) < d(p H2 ,Pi), then i e H[, else i G H 2 , (12) 

where 

d(Pj,Pk) = AR(pj,p k ) = ^ A(()(p j ,p k ) 2 + Ar](pj,p k ) 2 . (14) 

A(p(pj,Pk) and Ai](pj,pk) are the azimuthal angle difference and the pseudo rapidity 
difference between j and k, respectively. Finally, we reassign Hi = {H[,fi} and 
H 2 = {H 2 ,f 2 }, and repeat the second step, i.e., the assignment of jets and leptons. 
The final step continues until the assignments to Hj (j = 1,2) converge. 

It would be useful to combine a hemisphere analysis to the mass measurement 
of qi- The biggest benefit is that the combinatorial backgrounds are reduced by per- 
forming the analysis on inside of each hemisphere. As mentioned above, since multiple 
jets and also r-jets are expected in a SUSY event, a reduction of the combinatorial 
backgrounds, keeping the number of the signal as many as possible, will be effective. 

Next we discuss the left-handed squark mass measurement by using another decay 
mode of qL, that is, q~L q'W ± — > q'v T fi. This is interesting because this analysis 
provides information on the left-handed squark mass, which is independent of the 
analysis through qi — > qW° — > qrfi. The signal consists of one jet, one stau and also 
v T . The momentum of v T from W ± decay, denoted as p\^ \ is not directly measurable 
but needed for the mass measurement. For the momentum reconstruction, we consider 
the event 

Ml -+ (qx^q'W*) -+ (grf 1 )(gV r f 1 ), (15) 

where x° (<?) stands for either B (q R ) or W° (c[l)- We first apply the hemisphere anal- 
ysis and separates reconstructed objects into two hemispheres. The first hemisphere is 
required to contain one jet, one r-jet and one fi. The squark reconstruction through 
q — > qx° — > qrfi followed by r — > u T j r has been discussed above, and the kinematics 
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is fully known. From Eq. (JZJ), the momentum of v T from r decay is written as 

pW = (r-l) Pir . (16) 

Now we move on the reconstruction of the decay chain — > q'W ± — > q'p T fi. The 
second hemisphere is required to contain one jet and one f\ but no r-jets. With the 
assumption that the observed missing transverse energy is equal to {pi} + p^} ^)t, 
we obtain x- and y-components of p^} \ Furthermore z-component of pi}} ^ is deter- 
mined from the condition of 



m% ± = ( P f ±} + P,J 2 , (17) 

with a two-folding ambiguity. In this way, the reconstruction of p}} ^ is possible, but, 
it would be useful for a background reduction to impose the condition 



M T (pf ±) ,^ 1 )<m 



W ± i 



where Mt is defined by Mt(p,p') = ym? + m' 2 + 2(EtE! t — p T ■ p' T ). Then, we 
would gain information on the q^ mass from the distribution of 



m^ l = vfe+pir ] +p^) 2 - as) 

Now, we introduce a sample MSSM-spectrum and perform a Monte Carlo analysis 
to demonstrate how the method described above will work. The underlying model is 
assumed to be the minimal gauge-mediated model with parameters: 

A = 60TeV, M mess = 900 TeV, iV 5 = 3, tan/3 = 35, sign(/i) = +, (20) 

where A is the the ratio of the F- and A-components of the SUSY breaking field, 
M mess is the mass scale of messenger fields, A^ 5 is the number of messenger multiplets 
in units of 5 + 5 representation, tan is the ratio of the vacuum expectation values of 
two Higgs bosons, and sign(/x) is the sign of the SUSY invariant Higgs mass parameter. 
Low-energy mass spectrum is calculated by ISA JET 7.64 [29], in which we use the 
top-quark mass of 171.3 GeV. Superparticle and the lightest Higgs boson masses are 
shown in Table [TJ This parameter point has been studied in early work [SU] , and it 
has been shown that B, W° and q~R masses are well determined. Hence we assume 
that these masses have been measured. We further assume that the mass of the lighter 
chargino, W ± , has been measured. Branching fractions relevant to the following study 
are as follows: BR(g L qW°) ~ 23%, BR(u L -> dW + ) ~ 49%, BR(d £ -± uW~) ~ 
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Particle Mass [GeV] 
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1309.39 




1231.70 
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1183.97 


d L 


1234.28 


d R 


1180.19 


h 


1082.85 


k 


1195.08 


h 


1145.24 


h 


1185.83 


$1 


388.05 


h 


396.19 
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402.57 


v T 


383.80 


en 


193.39 
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148.83 




239.52 


xl 


425.92 


xl 


508.41 


x\ 


548.67 


xt 


425.45 


xt 


548.43 


h 


115.01 



Table 1: Masses of the superparticles and the lightest Higgs boson h. 
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while qn decays into qB with almost 100% branching fraction. For the lighter 
neutralinos, BR(B — > ^ff) ~ BR(iy° — >■ r ± f 1 T ) ~ 31%, and for the lighter chargino, 
BR{W + -> v T ft) ^ 61%. 

The SUSY event generation and the hadronization process are simulated by HER- 
WIG 6.510 [301 EI]. Events generated are passed through the fast detector simulation 
package, PGS 4 [32], with a slight modification to treat a stable stau. We assume 
that the momentum resolution of stau track is same as that of muon, and the energy 
deposition of stau to calorimeters is negligible. For the stau identification, we impose 



# 2 The soft Wino mass parameter and /^-parameter are somewhat close at this parameter point. 
This causes the left-handed squark decay to the heavier chargino with a considerable branching 
fraction. 
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two conditions: (1) relatively small velocity, 0.4c < v < 0.91c; (2) large transverse 
momentum, px > 50 GeV. Such staus are identified as stau, while muons are distin- 
guished from stau under these conditions. A stau with large velocity, v > 0.91c, is not 
identified as stau but muon. The total SUSY cross section is 669.6 fb in pp collision 
with the center of mass energy of 14 TeV, and we generate all SUSY processes. We do 
not consider any standard model processes as background sources; In the long-lived 
stau scenario, SUSY events are characterized by the exotic slow-moving charged track 
of the stau. 

We first discuss the left-handed squark mass measurement in the sequential decay 
process, q^ —> qW° — > qrf\. Each of SUSY events should contain two stau tracks in 
the case of the long-lived stau. However, we require at least one stau in the signal 
event because f\ with large velocity is expected to be identified as muon. Notice that 
such a muon, which is a mis-tagged stau in reality, may have large px- Therefore, we 
regard the highest px muon-like track as the second stau track if there is only one slow 
stau in the event. This procedure ensures that we have always two stau candidates in 
a SUSY event. Then we perform the hemisphere analysis described above. We impose 
two conditions so that the hemisphere method works fine: 

IA) AR(p fl ,pr 2 ) > 2.0, where fi and f 2 are the first and the second stau in the event. 

IB) (p Hj ) T > 400 GeV for j = 1, 2. 

We regard objects in each of two hemispheres as daughters of a heavy superparticle 
produced in hard-collision. So, for each of the hemispheres, the following conditions 
are imposed to find out the signal: 

la) Exactly one stau. 

lb) At least one r-tagged jet with p T > 30 GeV and the charge opposite to the stau. 
lc) At least one jet with px > 200 GeV and m jet < 80 GeV. 
Id) No leptons with p T > 20 GeV. 

Now, we reconstruct the left-handed squark decay chain. The momentum of r is 
reconstructed from that of r-jet with the rescaling factor r^ (> 1), using W° mass 
information. In this analysis, a pair of r and fi from B decay is one of backgrounds. 
From Eq. ffTUj) and Eq. (fTTj) . we see that the relation of r^ > r§ always holds for any 
pair of {j T ,Ti) in the case of > rrig. Therefore, r-jet and fi from B decay with 
rg (> 1) will be reconstructed as if they are produced in W° decay because the signal 
condition of r^o > 1 is always satisfied due to the relation of r^o > rg. In order to 
reduce this background we impose not only r^o > 1 but rg < 1. Then we study the 
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Figure 1: Distribution of the invariant mass Mg L = (pj + r^ p T + pj 



distribution of the invariant mass 



(21) 



where pj (pj T ) is the momentum of the highest pt jet (r-jet). In Fig. (TJ we show the 
distribution of Mg L for the integrated luminosity C = 100 fb _1 . One can see the peak 
structure in the distribution, whose position corresponds to the squark mass. In order 
to determine the squark mass, we take the fitting function, 



f{M qL ) = Nexp[-(M, L - M^ cak) ) 2 /2a 2 ] + A{M, L - 



1L 



+ B, 



(22) 



where N, M^ eak ' ) , cr, A and B are parameters. Then, we obtain the value of M ( j^ eak ' ) = 
1225 ± 3.2 GeV while the underlying value is = 1232 GeV. The best-fit value is 
smaller than the input value. This may be because, for example, some of jet energies 
are leaked from the cone used in the jet reconstruction. Hence, such systematic errors 
could be corrected when the jet energy is appropriately calibrated. 

We also study the right-handed squark mass measurement, although it has already 
been discussed in early work, for a comparison. The conditions imposed in this analysis 
are same in the above. In Fig. [2], we show the invariant mass distribution of Mq R with 



rg > 1, where 



M; 



1R 



{Pj +^bPt +Pf 



(23) 
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Figure 2: Distribution of the invariant mass Mq R = (pj + r gp T + pf x ) 2 . 

We perform the fitting of the distribution by the fitting function introduced in Eq. fl22|) . 
and obtain = 1176 ± 2.3 GeV. (The underlying value is m 9 - fl = 1184 GeV.) 

Again, the small discrepancy between rriq R and M~^ eak ^ are seen, but it could be 
corrected if we take some sources of systematic errors into consideration. Comparing 
two figures, we see a clear difference between M^ eak ' ) and M~^ eak ^ . This is an important 
evidence that there really exist two types of squark, that is, one is the right-handed 
squark and the other is the left-handed squark with different masses. Also, notice 
that we see from the figures the number of the q~L signal is smaller than that of q~R 
signal. This indicates that the squark reconstructed with r^o is really the qi, because 
the smallness of the number of signal may be represented a small branching fraction 
of (Jl qW°. Inversely the squark reconstructed by using is implied to be qR. 
Of course, in order to confirm the discussions here, we should seriously consider the 
signal efficiencies in both analyses. 

There exists another decay mode of q~L into q'W ± , and this decay mode can be used 
for the mass measurement. We now discuss the left-handed squark mass measurement 
in the sequential decay process, qi —> q'W^ 1 —> q'v T fi. The signal consists of one jet, 
one stau and also v r from W ± decay. In order to obtain the momentum of this u T , as 
discussed above, we consider the event qqi — > ( ( ?X )( ( ?'^ ± ) ~~ >* (QT^il'^T^i)- We first 
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use the hemisphere analysis. One of hemispheres is required to contain: 
2a) Exactly one stau. 

2b) At least one r-tagged jet with pt > 30 GeV and the charge opposite to the stau. 
2c) Al least one jet with p T > 200 GeV and m jet < 80 GeV. 

If these conditions are satisfied, we calculate and r^o from Eq. ffTUl) and Eq. ffTTT) 
for each pair of (j, j T , fi). For finding out a squark candidate, we further impose: 

2d) If r B > 1, \Mg R - M^ eak) | < 40 GeV. 

2e) If r B < 1 and r Wo > 1, |M 9 ~ L - ild peak) | < 30 GeV. 



Here M 5l and M~ qR are defined in Eq. (f2T|) and Eq.(J23]), and we use M^ eak) = 1176 GeV 
and M~^ eak ^ = 1225 GeV, taking into account their systematic errors. We calculate 
the momentum of v T from r decay by using r§ if the former condition, 2d), is sat- 
isfied, or by using r^o if the later condition, 2e), is satisfied. Combining this recon- 
structed v T momentum with the observed missing transverse energy, we obtain x- and 
y-components of the momentum of v T from W ± decay. Now we are at the stage to re- 
construct the signal process, q~i — > q'W ± — > q'v T f]_, by analyzing the other hemisphere. 
We impose the conditions: 

2a)' Exactly one stau. 

2b)' At least one jet with p T > 200 GeV and m jct < 80 GeV. 
2c)' No r-tagged jets with p T > 30 GeV. 
2d)' No leptons with p T > 20 GeV. 
2e)' {pt ±] )T > 50 GeV. 
2f)' M T (p^ ±} ,p fl ) <m w± . 

We reconstruct z-component of p^f ^ by Eq. (|T7|) . Then, the squark mass is deter- 
mined by the peak analysis. As mentioned above, the two-folding ambiguity exist 
in this analysis, and we choose the solution which has a larger invariant mass 



#3 



111 



Fig. [31 we show the invariant mass distribution of Mq L = yj (pj + ^ + PfJ 2 , where 
we use the highest px jet in the hemisphere. The integrated luminosity is taken to 
be 300 fb _1 in this study. Using the fitting function given in Eq. fl22|) . we obtain 



# 3 If both solutions are included in the analysis, one will see a fake peak in the distribution, which 
comes from the wrong solution in the reconstruction of pu^ . In the present parameter point, a 
position of the fake peak is around Mq L ~ 1000 GeV. Such a fake peak will have been confirmed not 
to be a signal in reality, when we perform the Monte Carlo analysis using the M^ cak) ~ 1230 GeV 
and other measured superparticle masses as input parameters. 
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Figure 3: Distribution of the invariant mass Mq L = y (pj + pl^ ^ +Pn) 2 - 
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the values of M^^' determined by two independent measurements agree with each 



1232 GeV.) Interestingly, 
rrements agree with each 
other within the range of the statistical errors. On the other hand, we see a signifi- 
cant difference between Mg L = v/ (Pj + p^ ^ + Pn) 2 and Mg R = a/ (pj + r^p T + p^) 2 . 
This is another important indication that the squark reconstructed in the present 
analysis is q^. However, these results are obtained based on the simple detector sim- 
ulation, and also, they are sensitive to the choice of the fitting function, the fitting 
range and so on. Thus, a detailed study is needed to confirm the discussions. 

Mass measurements presented in this paper provide several information on the 
MSSM. Their masses themselves reflect SUSY breaking mechanism, such as the mes- 
senger scale or the SUSY breaking scale. After the mass measurements one could 
test the several SUSY breaking models by extrapolating the measured masses to high- 
energy scale using the renormalization group equations. Also, superparticle proper- 
ties could be determined by these analysis. For instance, q~L has two decay modes, 
qL — > q'W ± and — > qW°, and both of them can be used for the mass measurement. 
If we will have confirmed that the measured masses in two analyses agrees with each 
other, this is an evidence that the squark is a really left-handed one. Once squark 
chiralities have been determined, we can further investigate the other superparticle 
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properties by studying a variable which reflects superparticle natures; One of such vari- 
ables is the invariant mass of (q, r) system in the squark decay qu — > qB — > qr^r^ [21], 
and this variable is useful for the determination of the chiralities and spins of super- 
particles. 

In summary, we have studied left-handed squark mass measurements in the long- 
lived stau scenario at the LHC experiment. q~i has, generally, two decay modes and 
hence the number of the signal becomes small comparing to that of q~R. This is 
one of difficulties in q^ mass measurement. Another difficulty is a combinatorial 
background associated to the multiple jets and also r-jets, which is expected in the 
LHC experiment. In order to reduce backgrounds and find out signals, we introduce 
a hemisphere analysis in which we take into account the fact that there are two 
stau tracks in each SUSY event if the decay length of fi is enough long. We also 
propose some techniques to reconstruct the momentum of v T which is associated to 
q~L decay. Combining these materials, we have shown that m$ L is determined by peak 
analysis in both of two decay modes of the left-handed squark, q~i — > q\V° — > qrfi and 
q L -> q'W ± -> <?V r fi, at the LHC. 

Acknowledgments: The author would like to thank T. Moroi for fruitful discussions 
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